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Q ' Abstract 

in '. 

, We propose a new mechanism to explain the radio Arc and threads in the Galac- 

' tic Center by current sheets produced by local magnetic shears due to interac- 

Q^' tion of a moving cloud and the vertical field based on three-dimensional magneto- 



hydrodynamical simulations. Magnetic reconnetion and acceleration of cosmic-ray 



|H I electrons in the current sheet will result in high contrast of radio emissivity inside 

^ ■ and outside the Arc and threads. 

■ 
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1 Introduction 

The Galactic Center is full of spectacular radio features, among which the non- 
thermal radio Arc is the most prominent unusual structure, believed to be highly 
ordered magnetic field vertical to the galactic plane. Each magnetic tube is as thin 
as one parsec or less and as long as a few tens of parsecs (Yusef-Zadeh, Morris, M., 
Chance 1984; Yusef-Zadeh, Morris 1987; Morris, Yusef-Zadeh 1985). Besides the 
Arc, there are numerous vertical radio filaments, called radio threads, also believed 
to be of magnetic origin (Anantharamaiah, et al. 1991; Lang, et al. 1999; LaRosa, 
et al. 2000; LaRosa, et al. 2004). 

Magnetic fields, as strong as 0.01 — 1 mC, have been inferred for the radio Arc 
and threads from their large Faraday rotation measures (Lang, et al.l999; Tsuboi, 
et al.l986; Sofue, et al. 1987; Yusef-Zadeh, et al. 1997). Their straight structures 
indicate that they are composed of highly ordered magnetic fields perpendicular 
to the galactic plane. Although large-scale vertical magnetic fields appear to be 
common in the central regions of spiral galaxies (Tsuboi, et al. 1986; Sofue, et al. 
1987; Sofue, Fujimoto 1987), it has been a mystery why the Arc and threads are 
shining so locally, exhibiting high contrast of brightness to the ambient regions. 

A couple of models have been proposed for their origin. In a spark model for the 
radio Arc (Sofue, Fujimoto 1987), ejection of magnetized gas from Sgr A, observed 



as the thermal filaments, hit the vertical field in galactic rotation, and reconnection 
of the fields in the Arc and thermal filaments results in acceleration of cosmic- 
ray electrons to radiate nonthermal radio emission. On the other hand, a galactic 
wind model (Dahrlburg, ct al. 2002) suggests that interaction of the wind from the 
nucleus with a gas cloud causes wake instability, which grows to produce filamentary 
structures perpendicular to the galactic plane, mimicking the radio threads. The 
former theory, however, cannot explain the number of threads. The latter theory 
may explain the morphology of the threads, but cannot be applied to the strong 
Arc feature and is also difficult to eplain the high contrast of radio brightness in the 
threads. 

In this letter, we propose a new mechanism to explain the Arc and threads by 
current sheets produced by local magnetic shears due to ineteraction of a moving 
cloud and the vertical field. Magnetic reconnetion and acceleration of cosmic-ray 
electrons in the current sheet will result in high contrast of radio emissivity inside 
and outside the Arc and threads. 

2 Three-Dimensional Non-axisymmetric MHD 
Simulation of Cloud- Magnetic Field Interaction 

Vertical magnetic fields are often observed in spiral galaxies, including the Milky 
Way, and will be a common structure in the central regions (Sofue, Fujijmoto 1987). 
If there exists an accreting gas disk, the vertical field lines are twisted, and form a 
large-scale magnetic jet, and various simulations have been made in axisymmetric 
scheme (Uchida, et al. 1985; Uchida, Shibata 1985; Shibata, Uchida, Y. 1986, 1987). 

On the other hand, the distribution of molecular gas in the Galactic Center is far 
from axisymmetric, but is significantly shifted to positive galactic longitudes, and 
is very clumpy (Bally, et al. 1987; Oka, et al. 1998). Longitude- velocity diagrams 
also show non-axisymmetric kinematics, which is partly induced by accretion of gas 
by a bar structure (Sawada et al. 2004) . Such non-axisymmetric motions of gas 
will cause significant local disturbances of the magnetic fields, which may only be 
treated by a three-dimensional MHD calculations. 

In order to examine how a locally moving gas cloud disturbs the ordered magnetic 
field, we performed three-dimensional MHD simulations. The calculations have been 
carried out using the CIP-MOC-CT scheme (Evans, Hawley, 1988; Yabe, Aoki, 1991; 
Yabe, et al. 1991; Stone, Norman 1992; Kudoh, et al. 1999; Kigure, et al. 2004). 
The CIP scheme can accurately follow the contact discontinuity, or the surface of 
the cloud in the present case. This scheme has been developed to three-dimensional 
cylindrical code, and is now applicable to local disturbances of field lines near the 
origin of the coordinates (Kigure, H., and Shibata, K. 2005 in preparation). The 
numbers of meshes are 171 x 100 x 480 in the radius, azimuth, and height (r, (f), z) 
directions. 

We assume that the ambient magnetic field is at rest with respect to the rest 
coordinate frame, and the filed lines are perpendicular to the galactic plane defined 
by (r, 2: = 0). The gravitational force was added by putting a point mass M at 
the Galactic Center coinciding with the origin of the coordinates, which mimics the 
central massive core (Takamiya, Sofue 2000) including the black hole (Genzel, et 
al. 2003; Ghez, et al. 2005). We put a gaseous corona around the central mass at 



hydrostatic equilibrium, but neglect its rotation. 

The initial field strength is taken to be uniform at B = Bq. A gas cloud with a 
radius is put at an initial position of (r, (f),z) = (1, 0, 0). The cloud is given of a 
Keplcrian initial velocity {vr,V(i),Vz) = (0, 1,0). The gas cloud is assumed to have 
an initial density distribution represented by 

P = Po[l + cos(7r|x-xo|/rw)]/2 (1) 

if |x — xqI < with xq = {r,4>,z) = (1,0,0), and we took = 0.3. The density 
distribution of the corona around the center is given by 

p = /,cexpHro/(r2 + z2)i/2_i|]^ (2) 

where ro is the unit length. The parameter a is defined as a = tV^q/^c' where Sc 
is the sound velocity in the corona, Vko = {GM /roY^"^ is the Keplerian velocity at 
radius r = tq. Here, pc is the coronal density at radius ro, a = 1.0 and pd po = 10~^, 
where po is the initial cloud density at the cloud center {r,4>,z) = (1, 0, 0). 

The plasma-/?, the ratio of gas pressure to magnetic pressure defined by ^ = 
2s^/7V^, is taken to be 0.4 in the cloud center, and is of the order of ~ 10~^ in the 
ambient coronal region, where 7 = 5/3 is the specific heat ratio of the gas. Here, s 
is the sound velocity at the clod center. We assume a free boundary condition. The 
gravitational and kinetic energy densities in the cloud are much greater than the 
magnetic field energy density. Hence, the field lines are nearly passive to the motion 
and distortion of the cloud. Initially we take u^/u^ ^ Hm/itk = (Klo/^o)^ = 
5 X 10~^, where Va.o and Vko are the Alfven velocity and the velocity of cloud. 
Here, Ug = GMop/r, = pv^ and Ura = B"^ [Sir. On the other hand, we 
have (^Ac/Mio)^ 5 for the coronal gas, and therefore, the magnetic field is not 
strongly disturbed by the corona. The gas pressures between the corona and cloud 
are assumed to be balanced, so that (sc/s)^ = Po/Pc ~ 10^; where Sc is the sound 
velocity in the corona. 

Figure 1 shows the result of the three-dimensional MHD simulation. Vertical 
lines are magnetic lines of force. The (x, z) plane is indicated by a grey-colored 
panel, with the grey brightness being proportional to the coronal gas density. A 
gas cloud is rotating around the center in Keplerian motion. The top panels show 
the time evolution of field lines from t = to 2, and the bottom-left panel is at 
i = 3. Here, time t is normalized by P/2'k with P being the rotation period of the 
cloud around the center. As the cloud moves, the magnetic field lines are wound, 
and twisted locally due to the cloud motion as well as the gravitational distortion of 
the cloud shape. Note that the viewing point of figure 1 is rotating with the cloud 
motion, so that the (x, z) plane is apparently rotating in the figure. 

— Fig. 1 — 

According to the twisting distortion of the field lines, a current sheet is produced, 

along which magnetic rcconnection is supposed to take place. The bottom-left panel 
of figure 1 shows an equal-currcnt-dcnsity surface with J = 1.5 at t = 3. Here, J 
is the current density calculated by J = roti?, and is normalized by Jq = bo/vQ = 
[poV^f^^f^ /ro, where 60 is the non-dimensionalized initial magnetic strength. We 
emphasize that the thus formed local enhancement of current density produces a 
feature mimicking the radio Arc and thread. The bottom right panel shows the 
radio Arc as reproduced from VLA observations at 21 cm (Yusef-Zadch, Morris 
1987). 



3 Magnetic-Reconnection and Current Sheet 
Model for Arc and Threads 



Based on the simulation, we propose a new model to explain the origin of radio Arc 
and threads in the Galactic Center, which we call the magnetic-reconnection and 
current-sheet (MRC) model. 

Suppose a gas cloud with different velocity from that of the field lines hit the 
vertical magnetic lines of force, such as due to infalling cloud, ejection from the 
center, or shock waves. The cloud will locally twist the vertical magnetic field. 
Figure 2 illustrates how a moving cloud influences the magnetic field lines. The 
locally twisted field lines, then, produce a magnetic shear between the twisted and 
the ambient fields. This magnetic shear produces a current sheet, along which 
magnetic reconnection would occur. This mechanism is similar to the nano-flare 
model for solar-coronal heating (Parker 1988). In the present case, the twisted field 
looks more like a bunch of lines of force than a sheet, and, hence, the current sheet 
would be more like a "current thread" , as figures 1 and 2 show. 

— Fig. 2 — 

The fiux density of the released energy per unit volume along the current sheet 
can be expressed by the Poynting flux, 

/ ~ {B,B^/A'k)V 

where i?^, i?^, and V are the magnetic strengths in the vertical and azimuthal 
directions, and the velocity of gas, respectively. 

The vertical component is on the order of the initial value, B^ ^ B^q. On 
the other hand, the initial value of is almost zero, while it is created during the 
twist to attain a finite value i?^, which wc denote as B^ ~ C-^zO- Hence, the flux of 
released energy along the current thread is of the order of 

/ ~ ^(S,oV47r)y. 

In the present numerical calculation, we have ^ ~ 0.1. This flux should be compared 
with the one in the ambient region, 

/o ~ (B,o-B<^o/47r)Vb. 

Thus, the ratio of the released fluxes inside and outside the current thread is of the 
order of 

r/~///o~e[5W50o][W]- 
The ratio rj will attain an extremely large value, because the ambient region has 
almost negligible value of B^^, as there exists neither a magnetic twist nor a current 
sheet. 

The energy released in the current sheet will be used to accelerate cosmic-ray 
electrons, which interact with the magnetic field and radiate synchrotron emission. 
We emphasize that the radio emissivity in the current sheet (thread) increases sig- 
nificantly even without large amplification of the field strength because of the large 
•q value. Let the fraction 5 of the Poynting flux / be transformed to the synchrotron 
emission. Then, we obtain the radio flux of 



The surface brightness of radio emission at frequency v is, then, estimated to be of 
the order of 

T. = f/u^ 5iBlV/AiTv. 

If we take Bq ~ lO^/itG, V ~ 100 kms s~^, and ~ 1 GHz, we obtain S ~ 
5^10~^^ erg s~^ cm~^ Hz~^, or S ~ (5.^10"^^ W m~^ Hz~^. This leads to brightness 
temperature at 1 GHz of T;, = 5^A^S/2/c 10^5^ K at around i/ ~ 1 GHz, where 
\ = c/u IS the wavelength and k is the Boltzmann constant. 

The radio brightness for the radio Arc and threads has been observed to be of 
the order of ~ 0.1 K at ~ 1 GHz (Yusef-Zadeh, Morris 1987; LaRosa et al. 2004). 
Here, ^ is approximately '^0.1 from the simulation. The transformation efficiency 
5 of Poynting flux to radio emission is a subject to sophisticated treatments of 
magnetic reconnection and acceleration of cosmic rays. Instead, we, here, estimate 
a possible value of the paprameter 5. In order for the formulated radio brightness 
to agree with the observed brightness, we require only a very small efficiency on the 
order of 5 ~ 10~^. 

We have discussed the MRC model for one particular cloud. However, due to 
barred galactic shocks as well as angular momentum loss by large-scale magnetic 
field twist, a number of clouds are infalling to the Galactic Center. Recurrent 
interaction of gas clouds and clumps will, therefore, result in many "threads" in 
the Galactic Center region, as indeed observed (LaRosa et al. 2000, 2004). Among 
the many threads, the strongest interaction between a cloud and magnetic field is 
presently observed at the radio Arc. In fact, the field lines in the Arc are interacting 
with a dense gas cloud (Tsuboi, et al. 1997). Interaction with the thermal filament 
extending from Sgr A with the Arc (Sofue, Pujimoto 1987) may also infer local 
twist. Siich disturbances will be caused not only by accreting gas clouds, but also 
by hydrodynamical waves and shocks triggered by various activities in the Galactic 
Center region, such as the star formation, gas ejection from the center, etc.. Analyses 
of the shapes, distribution, and frequency of radio threads may give information 
about hydrodynamical conditions in the Galactic Center region. 

Numerical computations were carried out on VPP5000 at the Astronomical 
Data Analysis Center of the National Astronomical Observatory, Japan (project 
ID: rhk05b and whkOSb), which is an interuniversity research institute operated by 
the Ministry of Education, Culture, Sports, Science, and Technology. 
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Figure Captions 



Fig. 1. Three-dimensional, non-axisymmetric MHD simulation of interaction of a 
gas cloud with vertical magnetic field in the Galactic Center region. The field is 
embedded in a hydrostatic gaseous halo around a central mass. Upper panels show 
the evolution of magnetic lines of force at i = 0, 1 and 3. The moving cloud is 
disturbed by the gravitation, and locally twists the magnetic field. For the units 
of time and linear scale, see the text. The lower-left panel shows an equal-current- 
density surface of the current sheet (current thread) at J = 1.5 for t = 3. The 
lower-right panel is the radio Arc (vertical structures) near the Galactic Center (the 
brightest spot) (Yusef-Zadeh, Morris 1987). 

Fig. 2. Schematic illustration of a scenario for the origin of Galactic Center radio 
threads, due to local twist of vertical magnetic field lines. Reconnection of the 
field lines along current sheets (threads) accelerates cosmic-ray electrons to emit 
synchrotron radiation. Inserted in the bottom-right corner is the observation of 
90-cm radio threads (LaRosa et al. 2000). 
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